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Abstract. The idealized theory of quantum vacuum energy density is a beau- 
tiful application of the spectral theory of differential operators with boundary 
conditions, but its conclusions are physically unacceptable. A more plausible 
model of a reflecting boundary that stays within linear spectral theory confines 
the waves by a steeply rising potential function, which can be taken as a power 
of one coordinate, z a . We report investigations of this model with consider- 
able student involvement. An exact analytical solution with some numerics 
for a = 1 and an asymptotic (semiclassical) analysis of a related problem for 
a = 2 are presented. 



1. Introduction 

The Casimir effect [9l [8j 129L I27L [6] is an observable attraction between neutral 
electrical conductors. Its mathematical charm is that, at least for perfect conduc- 
tors, it can be attributed to and calculated from the geometrical dependence of the 
energy of the quantized electromagnetic field in the region between the conductors. 
As often in quantum field theory, naive calculations give an infinite answer. Sub- 
traction of the zero-point energy of each field mode renders the local energy density 
finite but leaves a nonintegrable singularity at the boundaries. This divergence, 
because it can be regarded as permanently attached to the conductors, does not 
interfere with the calculation of the force of attraction between rigid bodies. How- 
ever, a more accurate representation of the physics near the boundary is needed 
[3l 120] to understand situations where the boundary can deform, such as the cel- 
ebrated case of an expandable sphere (7J. Furthermore, the energy density (more 
completely, the stress tensor) of the field should act as a source of the gravitational 
field in general relativity, so localized infinities within it are not physically plausible 
|10j . For more detailed background information see |14l 1151 llll 1171 121] . 

It is universally agreed among physicists that the root of this problem is that the 
idealization of a "perfect conductor" is inapplicable to very-high-frequency modes of 
the quantized field. A full treatment of the physical problem, including modeling of 
the charged particles inside the conductors, would take us out of the spectral theory 
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of self-adjoint linear partial differential operators into difficult, nonlinear condensed- 
matter physics [4]. One might hope that an ad hoc cutoff of the contribution 
of high-frequency modes would preserve the qualitative essence of the physically 
correct solution, and indeed a simple exponential cutoff leads to both tractable 
calculations and physically plausible results |17l 12 lj , while placing the topic firmly 
within the study of the asymptotics of integral kernels (Green functions) associated 
with the operator concerned, in the grand tradition of spectral and geometrical 
asymptotics. (Most of the issues of principle in vacuum energy are adequately 
addressed by studying a scalar field with Dirichlet boundary conditions instead 
of an electromagnetic field with perfect-conductor boundary conditions, and our 
discussion has now lapsed into that setting.) 

Unfortunately, a close examination of the stress tensor (in particular, energy 
density and pressure) predicted by the theory of [IT] reveals that the expected 
energy-balance equation (sometimes called principle of virtual work [5]V 



describing the change in total energy when a boundary acted upon by pressure p is 
moved a distance h, is not satisfied |16j . The root of this problem appears to be that 
the degree to which a particular normal mode is affected by the cutoff depends on h, 
so that derivatives of the cutoff function contaminate fundamental relations like 
(jl.f p . The exponential cutoff on frequency amounts, after analytic continuation, to 
a time-direction separation of the space-time points that are the argument variables 
of the integral kernel for the wave equation of the problem. Physically plausible 
results can be achieved for the various tensor components in various scenarios by 
choosing a specific space-time direction for the point separation in each case, but 
such an ad hoc procedure cannot be regarded as logically satisfactory for the long 
term. 

Our present goal is to replace the reflecting boundary and the cutoff with 
a smooth potential rising to infinity. The potential can be thought of as a static 
configuration of another scalar field (presumably of very high mass) . After some well 
understood renormalizations in the bulk, this model should yield a finite vacuum 
stress tensor without any cutoff. As an internally consistent physical system without 
ad hoc modification, it is expected to satisfy the physically required energy-balance 
relation. Yet in the limit of a very steep wall it should approach in some sense the 
vacuum stress of the problem with Dirichlet boundary. If its predictions resemble 
those for hard-wall calculations with certain point separations, those (technically 
easier) ad hoc methods will be vindicated and can be used with confidence in other 
situations. 

Here we report progress on this problem achieved during the spring semester 
of 2010, while J. Wagner held a Visiting Assistant Professorship in Mathemat- 
ics at Texas A&M University. He worked closely with Professor S. Fulling and 
four research assistants supported by his NSF grant: mathematics M.S. students 
J. Bouas and F. Mera and undergraduates (physics-mathematics double majors) 
C. Trendafilova and K. Thapa. This brief period was notable for intensity of col- 
laboration. 



(1.1) 
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Figure 1. Graph of the potential v(z) for the two simplest choices 
of a. 



2. The model 

Consider the potential function (1 < a £ R; r = (x,y, z) e R 3 ) 

fO. z<0 
(2.1) v(r) = I ( \« 

Mi) , *>0. 

It is characterized by a and the single length scale 

A 

We ordinarily fix Ao = zq = 1 and let a vary (and suppress the arguments x and y) . 
Thus v(l) — 1 for all a and the potential forms an increasingly steep wall near 
z = lasa^oo (Fig. [TJ). 

The scalar field is an operator-valued distribution satisfying 

d{x ) 2 

Because the PDE (|2.2[) is linear and invariant under translation in time (a; ), it 
has a standard reduction to a classical eigenvalue problem: tp(x , x, y, z) is a linear 
superposition (an integral, since the spectrum in our model will be continuous) of 
normal modes of the form (j>j{r)e ±%u s x with operator coefficients. For the details 
we refer to [291 117j . The physical quantities of interest are (the vacuum expectation 
values of) the components of the energy-momentum tensor, especially the energy 
density 



(2-2) 757z£a=VV-»(*)v. 



(2.3) T 00 (r) = \ 



\dx 



ipW 2 ip + vtp 2 



These expectation values can be expressed as derivatives of the cylinder kernel 
(Poisson kernel) of the differential operator on the right-hand side of (|2.2[) , which 
can be defined by the PDE 

(2.4) (|L + V 2 - v(z)) T{t, x, y, z, z') = 2S(t)S(x)5(y)5(z - z') 



dt 2 

or by the eigenfunction expansion 



(2.5) T(t,r,v') = - J dM)u3f x il> j {T)<l> i <j>)e- u '\ 

^We take h = 1 = c, curvature coupling constant £ = j, and metric signature goo < 0. 
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where fi is a properly normalized spectral measure over the index set. (T is an 
analytic continuation via t = i{x — x') of the Green function of the wave equation 
(|2.2j) . In (|2.5[) we take advantage of the translational invariances of (|2.2I) to set 
x' = y' = without loss of generality.) Now (formally) 



1 fl 2 T 



t— 0,a:— 0— y,z'- 



and there are similar formulas (involving derivatives with respect to the spatial 
variables) for the three pressures, P z — T zz etc. As usual, we must refer to the 
references for complete explanations. 

In (|2.6[) and similar formulas, the integral kernel is being evaluated at equal 
arguments ( "on diagonal" or "in the concidence limit", depending on whether one 
is in mathematics or physics literature). In studying a new model it is common to 
look first at the diagonal value of T itself (without any derivatives) , which has the 
physical interpretation of vacuum expectation value of the square of the field (p. 



3. Analytical solution 

3.1. Eigenfunctions. Solving (|2.2p by separation of variables is a standard 
exercise in quantum scattering theory. The formal index j in (|2.5[) is a triple 
(k x ,k y ,p) with = (k x .k y ) 6 R 2 and p 6 (0,oo). The frequency cjj is the 
positive solution of w 2 = kj_ 2 + p 2 . The eigenfunction is 

(3.1) &(r) = (2^)- 1 e ik - r ^ p (z) 
where 

(3.2) J^ +w(z) _/) 0p ( z ) = o. 

When z < 0, it must take the form 

(3.3) 4> p {z) - sin[pz - S(p)] 

for some real phase shift S(p) (not to be confused with the Dirac S in (|2.4[) and 
(13.241) ). The normalization factors in (I3.1[) and (|3.3[) guarantee that the spectral 
measure /i is 3-dimensional Lebesgue measure. 
When z > 0, <j> is best expressed as 

(3.4) p (z) = C(p)P Q (z/z,(zp) 2 ) 
where the function P a (z,E) is a solution of 

(3.5) (^-^+z a -£)p Q (z,£) = 

that vanishes as z — > oo. For small, integer a the solutions are known as Airy 
functions and parabolic cylinder functions: 

(3.6) P 1 (z,E)<xM(z-E), P 2 (z,E)kDi (e _ 1) {V2z). 

For a hard wall at zo, we have P 00 (z,£ l ) ex sin[v / P(z — z o)]- Henceforth we take 
z = zq = 1 so that E — p 2 . 
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The solutions and their derivatives must match at z — 0. Thereby C and, more 
importantly, 8 are determined: 

(3.7) ^ m = ~P^y 

2 1 



(3 ' 8) C{P) it P a (0,P 2 ) 2 +p- 2 ^(0,p 2 ) 2 ■ 

Even in the cases (|3.6I) these formulas do not lend themselves to exact evaluation of 
the integrals for energy density and pressure, so further approximation or qualitative 
analysis is needed. 

3.2. Asymptotics. When p = the solution of (|3.5|) is known (a modified 
Bessel function): 



(3.9) 



P a{ z,0) = z^K^ (^^) 



For small p the solution can therefore be constructed as a perturbation expansion: 

(3.10) P a (z,E) = P a (z,0) + EP^(z) + --- . 

This process requires constructing the Green function (resolvent kernel) for the 
nonhomogeneous unperturbed equation and applying it iteratively. In this way we 
find that 

(3.11) w _ Ko + 2) *r(2±§)r(2±i)~W). 

Fortunately, the Bessel-function integral needed to find the 0(p 3 ) term can be 
evaluated in closed form, but we do not report the result here. 

At large p one can construct a WKB (semiclassical) approximation: 



(3.12) <t> P (z)~ [p 2 - v(z)}—± co 

where a = p 2 / a is the turning point. It follows that 



vV -v(z)dz- - 



S(p) 



/ \J p 1 — v(z) dz H mod ir. 

Jo 4 



Closer examination shows that the "mod 7r" can be ignored and the integral eval- 
uated as a beta function: 

(3.13) f(p) = i P i+»/-bC| ^)+7 + °W- 

In summary, we have for a — 1 (the Airy function) 

(3.14) S( P ) 



P 3 2 / 3 r(|)/r(|), P ^o, 



and for a = 2 (the parabolic cylinder function) 

'2 P r(|)/r(|), P ^o, 



(3.15) <S(p) 

I " 4 1 4 

We are interested in the power potential (|2.1[) only as a convenient model with 
suitable qualitative properties. Since the function 5(p) completely encodes the 
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influence of the potential on the field in the potential-free region, it is tempting to 
forget the potential and study the class of models parametrized by functions i5 in a 
suitable class. The asymptotic relations (|3.1ip and (|3.13[) give some idea of what 
an allowed S must look like, but otherwise the inverse problem of determining v, or 
even its basic qualitative properties, from a given S is wide open, as far as we know. 
Later we shall show evidence that slight changes in d can produce unacceptable 
results. 

3.3. The renormalized cylinder kernel. To exploit the symmetry between 
t and rj_ in this problem, we introduce another layer of Fourier transformation into 
(1231) : 

(3.16) T(t,r ± ,z,z') = —±— [ dudki_r dpe™e*^Mz)t{u,lL X ,p), 

(2tt)-V^ J r3 J q 

-2 <j> p (z') 



(3.17) T(u,\L±,p) = 



(2tt) 3 / 2 uj 2 + k ± 2 + p 2 



where u has now become an independent parameter. 

The integral over R 3 can be done by standard methods, resulting in 

— If 00 
(3.18) T(t,r ± ,z,z') = - — dpY{s,p)cj> p {z)ct> v {z% 

Z7T ' 







Y( 8 ,p) = —, s^^t 2 + \v ± \ 2 . 
s 

We shall concentrate for now on the potential- free region, z < 0, where 
— If 00 

(3.19) T = -—\ dpY(s,p)sm(pz- S(p))sm(pz' - S(p)). 

Upon converting the product of sines to a sum of cosines in the standard way, 
one sees that the first term is just the "free" kernel that would apply in R 4 if the 
potential were not there: 

^ 3 ' 2 °) T = ~2^t 2 + r ± 2 + (z- z') 2 

1 f°° 

+ 2^Jo * r ( s 'P) C0S (p( z + z ') - 26 (P)) 

— -^free "T" T rcn . 

For a hard (Dirichlet) wall at z — zq we have 8{p) = zqp and hence the well 
known image solution, 

1 1 



-Sp 



(3.21) T ren = 



2tt 2 t 2 + r ± 2 + (z + z'- 2z ) 2 



Before continuing it is instructive to take a close look at this case. Tf ree is, of 
course, singular on the diagonal (t = 0, = 0, z' = z) and only there. The 
singularity makes it impossible to pass to the diagonal directly in (|2.6|) . but that 
is also unnecessary: Tf roc is present in all problems, including empty space (where 
rpfiv j s na turally defined to be zero), and hence is physically meaningless. One 
expects to isolate and discard it before implementing (|2.6[) . which is applied only 
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to the "renormalized" kernel, T ron • The latter, as given in (|3.21[) . is nonsingular in 
the region of physical interest; on diagonal one gets (with our convention zo = 1) 

(3 - 22) Trcn = 8^ (z-l) 2 

and similar formulas (proportional to (z— 1)~ 4 ) for the energy density and pressures. 
In the present case f|3 . 19[) is applicable over the whole range — oo < z < 1, and 
(|3.22p gives the expectation value of the scalar field right up to the wall (where it 
develops a nonintegrable divergence). This much is totally standard and familiar 
to all workers in the field of vacuum energy. However, let us go back to the integral 
form of (EOT) contained in (pOO]) : 

— 1 f 00 e~ sp 
(3.23) T rcn = —j / dp cos(p(z + z') - 2p) , 

Z7T In S 



where s = yjt 2 + \r±\ 2 , and attempt to set t and r equal to before evaluating 
the integral. Although T rcn is well-defined except when z + z' = 2, in (|3.23[) we 
appear to have a double disaster: the denominator of the integrand is identically 
zero, and, moreover, even the integral of the numerator alone diverges because the 
exponential cutoff is lost. The resolution of this apparent paradox is that 

r °° 

(3.24) / cospz dp — tt8(z) 

Jo 

in the sense of distributions, and the Dirac distribution 8{z) is identically for z ^ 0. 
Thus the classically divergent numerator integral is equal to in the distributional 
sense in the limit s — > 0, and (|3.23[) is consistent with (|3 . 2 1 1) and (|3.22[) there. 
Returning to the general case, we are confronted by the integral 

— I f 00 e~ sp 

(3.25) T rcn = — / dp cos(p(z + z') - 2S(p)) . 

One might consider evaluating it numerically, given a trustworthy formula or ansatz 
for 5(p). However, in view of the previous special example, it is not surprising that 
the integral is poorly convergent when s = ^/t 2 + \r±\ 2 is small, which is precisely 
where we want it. In fact, we should be able to take s = and get a finite 
answer when z + z' > 0, but instead we have the same apparent infinities as in the 
Dirichlet case; and this time it is not obvious that the integral (without the factor 
s -1 ) vanishes distributionally to lowest order in s when z + z 1 > 0, although that 
must surely be true. 

To compound the problem, it appears that even this weak kind of convergence 
depends sensitively on S. Suppose that instead of (I3.23|) we had considered 

1 f°° e~ sp 
2t^ J d P—^- cos (p( z + z ')~ 2 P+§)> 

which one might naively think to correspond to the large-a limit of ()3.13j) (which 
actually is invalid for a — ¥ oo with fixed p) . This integral equals 

1 z + z' - 2 I 
" 2^2 V^ 2 + ri 2 t 2 + r ± 2 + (z + z> - 2) 2 ' 
Thus the cancellation that removes the divergence on the z axis in Q3.23P does not 
happen here. More generally, there is a genuine divergence for 5(p) = Ap + B 
unless B = 0. Naively one would think that the divergent boundary energy we are 
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studying is contributed by the modes of large p, and that therefore only the leading 
term in the WKB asymptotics (|3.13|) would be significant; the current example 
shows that that is not true. The fallacy in the reasoning is that high frequency can 
correspond to large at fixed p, as well as to large p. 

The sensitivity of (|3.23[) to a constant phase shift remains visible in the polar 
framework treated in the next subsection. The u integral in (|3.27[) or (|3.28j) in that 
case evaluates to a Bessel function J\ , which decays slowly as p — > oo; however, 
the outer integration yields (|3.22p by a standard formula found in handbooks and 
known to Mathematica. But the tiniest phase shift augments the Bessel function 
by a Struve function, which approaches a nonzero constant at infinity, so that the 
integral diverges unambiguously. 

3.4. Polar coordinates. Therefore, we recast the integration so that all high 
frequencies are treated on an equal footing. Any true divergence must come from the 
integral over high frequencies, since the eigenfunctions are smoooth and bounded. 

Abandoning the key formula (|3.25[) for now, we return to (|3.16|) . which in the 
notations Z = z + z' , s = (t, r±), v = (w, kj_), becomes 

(3.26) T rcn = — / dp dv^— - cos( P Z - 28(p)). 

4tt 4 J J R 3 v 2 + p 2 

(Note that s = |s| is the same s as before.) Now introduce polar coordinates in 
the space of variables (p, vi,V2,V3), with the main axis in the Z direction and the 
prime meridian through s (that is, S2 = = S3). After several steps of calculation 
one arrives at the new key formula 

(3.27) T len = — / dp du s~ x sin(spy 1 — u 2 ) cos^Zpu — 28(pu)). 

77 Jo Jo 

Because of the sine, the integrand is not singular, although it needs to be defined 
by a limit when s = 0. Ultimately we would like to take derivatives of (|3.27p and 
then pass to the diagonal, but for now we set s = and z = z' immediately: 

(3.28) T rcn (0,0,0,z,z) = — / dp dupyjl-u 2 cos(2z pu - 26 (pu)). 

71 Jo Jo 

(As previously remarked, this should yield the expectation value of ip 2 at z.) 

The change of variables has not eliminated the problem, but it has isolated 
it at the upper limit of a single improper integral. The convergence of (|3.28p is 
still delicate. Numerical and analytical investigations of it are ongoing. Were it 
not for the convergence issues, one could prove easily from (|3.28l) that the func- 
tion manifests approximately inverse-square decay resembling (|3.22| . but with the 
efffective wall position z — 1 replaced by z = c, where c is the coefficient of the 
linear term in p. lip . Numerical integrations have been performed in Mathematica 
for a — 1. Despite the instability of the highly oscillatory integrals, the results are 
qualitatively as expected, approaching (z — c)~ 2 /8tt 2 already for moderately large 
\z\ (Fig. [2]). 

As this paper was being drafted, K. Milton et al. reported (privately; cf. [28 ) 
an independent evaluation of T lcn for a — 1 using a different integral representation 
with better convergence properties. They find a weak (z~ l ) divergence in the energy 
density at the origin, which in hindsight is to be expected as diffraction from the 
coefficient singularity there (which will become less important for larger a). The 
corresponding singularity in T is of order zln \ z\, therefore not visible in Fig. [5] 
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FIGURE 2. (Solid) T rcn (|3.28[) as a function of z for a = 1. 
(Dashed) T rcn for a hard wall at c = 3 2 / 3 r(f )/r(§) w 1.37172. 

3.5. Conclusions. Our principal results are the Cartesian formula (|3.25[) . the 
polar formula p.27j) . and the latter's diagonal specialization ()3.28p . along with the 
formula (13. 7[) for 5 and the asymptotic analyses in the subsection following it. 
Formulas for energy density and pressure can be obtained by differentiating (I3.27|) 
and again passing to the diagonal. Granted the convergence of the integrals, it can 
be seen that T 00 and T 11 (the pressure in the x direction) are equal and opposite, 
so that (jl.ip is satisfied in the sideways motion of a plane partition perpendicular to 
the plane wall. This happy result, which fulfills the main motivation of the project, 
follows from the finiteness of the diagonal values (in the potential-free region) of 
the kernel T ren and its derivatives without the need of an artificial cutoff, since it 
is satisfied by the individual normal modes in the integrand. 

3.6. Inside the wall. Detailed calculations in the region z > have not yet 
been attempted, since we want to exploit fully the more elementary calculations in 
z < first. Also, one may reasonably consider the energy and stress in that region 
to be part of the wall, not part of the cavity containing the field. Nevertheless, 
in our setup the properly renormalized stress tensor inside the wall should still be 
finite and display the physically sensible behavior (jl.ip . (Most of the arguments of 
the previous subsection still apply here, or can be adapted.) 

On general grounds one knows that the renormalization will require subtraction 
of additional terms from T . In a nonsingular model the diagonal divergences in T 
and hence T^ v are determined [13] by the well known small-time divergences of 
the heat (or quantum) kernel [181 132) . Wherever the potential v(r) is not zero, the 
energy density calculated with the standard ultraviolet cutoff (t held nonzero) has 
the expansion 

(3.29) 7r 2 T 00 - f t- 4 - ±vt~ 2 + ±(v 2 - l^ 2 v) mi as t -> 0. 

The first term corresponds to the universal vacuum energy that has been removed 
by subtracting the free kernel at p.20j) . The other terms in (|3.29|) can be removed 
by subtracting higher-order terms in the small-i expansion of the exact T (obeying 
(|2.4j) with the potential included). Physically, such terms represent redefinitions of 
the constants in the equation of motion of the v field itself. 
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A technical issue that must ultimately be addressed here is the validity of 
()3.29p when v ^ C°° . In particular, for a — 1 (|3.29|) predicts a Dirac delta at z = 
from the term V 2 i>, and we will not trust the numerical coefficient until a direct 
calculation has been carried out. The issue here is the same one that arises in the 
heat kernel (and Gauss-Bonnet theorem) for the Laplacian in a region in the plane: 
The contributions of the corners of a polygon cannot be obtained by naively taking 
the limit of the effects of the curvature of a smooth boundary. This is an interesting 
question to be investigated in the future. 

4. Semiclassical analysis 

4.1. General remarks. Separation of variables, even when it is available, 
is not always the best way to obtain information about the integral kernels (and 
spectral measures) associated with a linear partial differential operator. Leading 
terms such as those in (|3.29[) are routinely found by direct construction of some 
kernel as an asymptotic series. Higher-order information can be obtained from 
terms in the quantum kernel (Green function of the time-dependent Schrodinger 
equation) corresponding to periodic orbits of the underlying classical mechanical 
system [TJ [2l 1191 I34j . (These terms are also present in the heat kernel but ex- 
ponentially suppressed. They produce oscillatory terms in the averaged eigenvalue 
density.) 

The construction of Green functions for the Laplace and Helmholtz equations 
in bounded domains in R™ (billiards) by reduction to integral equations on the 
boundary is well known. The counterpart construction for the heat equation is less 
familiar but available in the literature |22l 130] , What is seldom appreciated is that 
for the heat equation the solution of the boundary integral equation by iteration 
is convergent, because the integral operator has Volterra structure. One therefore 
has, in principle, an explicit construction of the solution. The Schrodinger equation 
has the same Volterra structure, so one expects again to have a convergent series 
solution. To implement this idea in a general context, Mera |25j has proved the 
following general Volterra theorem: 

Theorem 4.1. Let the kernel A(t,r) be (for each t and t in an interval I) 
a uniformly bounded linear operator A: B — > B, where B is a Banach space, and 
suppose that it has the Volterra property, A{t, r) = when r > t. Define the integral 
operator Q: L°°(I;B) -> L°°(I; B) by 

(4.1) Qm= I A{t,T)<t>{r)dT. 

Jo 

Then the Volterra integral equation 

(4.2) <f>-Q<P = f (f eL°°(I;B)) 

can be solved by successive approximations. That is, the Neumann series converges 
in the topology of L°°(I;B). 

The application of the theorem in any particular context reduces to showing 
that the operator family A(t, t) is uniformly bounded on a suitable space B. For 
Schrodinger equations this is a nontrivial task and requires supplementary technical 
assumptions. Here we are primarily interested in problems with potentials in R™. 
In that setting the key idea, due to Balian and Bloch [2 J , is to let the semiclassical or 
WKB approximation to the quantum kernel play the role played by the free kernel in 
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billiard problems, so that the role played by scattering off the boundary in billiards 
(or by scattering by the potential in standard time-dependent perturbation theory 
24, Ch. 7] 25J) is played here by scattering by a source that is essentially the 
residual error in the WKB approximation to the exact kernel. This construction is 
developed in 24 Ch. 8L 
The WKB kernel istf 

(4.3) G scl (x, t; y, 0) = (2™ft)-"/ 2 A e lS '\ 

where 



(4.4) 5(x,y,i) = / L(q(r),q(r))dr, L = 4q 2 — u(q) 

Jo 



is the classical action, a solution of the Hamilton-Jacobi equation, and the ampli- 
tude A is 



(4.5) A(x,y,t) = Jdet 



a 2 5 



If there is more than one classical trajectory q(r) starting at y at time and 
arriving at x at time t, the semiclassical approximation is a sum of such terms, 
possibly modified by Maslov phase factors (see next subsection) to keep track of 
places where the radicand in (|4.5[) has gone negative. 
Define a kernel Q by 

(4.6) Q(x,i;y,r) = ft 2 [A x A(x, t; y, r )] e ^(x,t;y,r) _ 
The corresponding operators Q and G sc i are related by 

(4.7) (-ihd t - h 2 V 2 + u(x))G sc i = 1 - Q ; 

that is, Q = 0(h 2 ) is the amount by which G sc i fails to solve the PDE for which it 
was devised. Thus, formally, 



or 

(4.9) 
where 



(4.8) G = G scl J2Q 3 , 

3=0 

G(x,t;y,r) = G scl (x, t; y, r) + f T(t, n)A(ri, r) dn + ■ ■ ■ , 

Jo 

(4.10) [r(t,r)^(r)](x) = / G scl (x,t; y, r)0(y, r) rfy, 

(4.11) [A(*,T)^(r)](x)= / Q(x,i;y,r)0(y,T)dy. 

Theorem 4.2. |24j Jn i/ie notation of the two foregoing paragraphs: Suppose 
that the following two hypotheses hold: 

(1) ||¥IIl-( W ») <00 - 

(2) f is o bounded operator from L 2 (R") to itself. 



2 In this section x and y are two different spatial points, not coordinates of the same point as 
earlier, and we reintroduce ti to make the structure of the semiclassical asymptotics clearer. To 
simplify the Schrodinger equation we take the mass m = i . 



12D. BOUAS, S. A. FULLING, F. D. MERA, K. THAPA, C. S. TRENDAFILOVA, AND J. WAGNER 




Figure 3. Paths with actions and amplitudes calculated in the 
text. Dashed: (5.13). Solid: (5.15) (heavy: (5.16-17)). Dotted: 
(5.19-20). 

Then the semiclassical operator A: L 2 (R n ) — > i 2 (R") is a bounded linear integral 
operator. It follows that the Volterra integral equation in the space L°° 2 (J; R") with 
the semiclassical kernel Q(yi,t;y,T) can be solved by successive approximations. 

This construction implements the Feynman path integral idea in a way different 
from the usual time-slicing approach. Each term in (|4.9|) is an integral over classical 
paths with j scatterings off an effective potential A A/ A. 

The determinant in (|4.5p is singular at caustics, where the mapping from initial 
velocity data (at y) to x ceases to be a diffcomorphism. One can expect both 
conditions (1) and (2) to be problematical if the orbit goes through a caustic, but 
we provide some evidence below that the situation is not as bad as one might expect. 
A way to go beyond caustics (if necessary) is provided by the Maslov theory 
as already implemented in a similar problem in 



4.2. The harmonic oscillator and the quadratic wall. It was natural to 
apply Theorem 14.21 to a power potential, with two motivations: to test the validity 
of the two hypotheses in the theorem in the context of a concrete problem, and to 
seek new information about the spectral density (and hence eventually the vacuum 
energy) for a soft wall. We have studied the case a = 2 in one dimension, 

^ fo if x < 0, 

( 4 - 12 ) v ( x ) = \ 1 2 2 - t Z n 

I jLU Z X Z if X > 0. 

The two transverse dimensions can be ignored because their contribution to the 
quantum kernel in dimension 3 is a trivial factor. 

Consider points x and y that are both in the potential-free region, and a time 
t > 0. There is always a force-free motion directly from y to x in time t. It is 
elementary to calculate for this direct path 

(4.13) S = ^ ~ V)2 , A 2 = ^, AA = 0, 
v ' At 1 2t ' 

so that (of course) the quantum kernel (|4.3|) is just that of a free particle. 

If Lot > 7r there is another classical path that enters the harmonic-oscillator 
region at 

(4.14) r = tl = _J_( t _I), 

x + y V uj J 

and reemerges after half a period, at to. = t\ + tt/lj. This path also contributes to 
the leading term, G sc \ , in (|4.9I) . (Bear in mind that such contributions are always 
of the schematic form Ae , which we shall not constantly repeat.) We shall see 
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that the action of the portion of the orbit inside the potential is 0, so the total 
action is just that of the two free motions at the ends: 

(4.15) S=f + A / = {X ( + V) \ , A 2 = - . 1 . , AA = 0. 

Note that the resulting term added to G sc i differs by a time translation (and a 
phase, since A 2 is negative) from the image term that would be produced by a 
hard wall. We shall show that the proper phase factor is —i. There is an apparent 
singularity in (14. 15)) at u)t = tt that deserves closer examination. 

Theorem 14.21 is formulated in [24] for a C°° potential. For (|2.1[) . in addition to 
(|4.15|) there are waves diffracted from the coefficient singularity at z = 0, but they 
become increasingly negligible with increasing a. 

Now consider x and y both inside the potential. From (|4.4I) and the relevant 
solution of the classical equation of motion one can reproduce well known formulas, 

(4.16) S(x, y, t) = 4g .^ ) [{x 2 + y 2 ) cos(^) - 2xy], 

(4.17) A 2 = — —. . , AA = 0. 
v ' 2 sin(wt) ' 

We need these formulas only for < ut < 7r; however, for the full harmonic 
oscillator potential on the whole real line it is well known |23|, 131) that the resulting 
(Mehler) formula for G sc i remains valid everywhere in space-time and gives the exact 
quantum kernel, with the caveat that (as suggested by the sign change in (|4.17jl ) the 
kernel must be multiplied by (— z) M where [i is the number of occasions when t has 
passed through an integer multiple of n/uj. (Conventionally one redefines A 2 and A 
to be always positive numbers and writes the Maslov phase factor (— i)^ separately.) 
At such a time there is a caustic; all the trajectories from y refocus at x = (— . 
But the kernel formula (|4.3j) reproduces there the original (t = 0) singularity, which 
is still a solution of the homogeneous Schrodinger equation. (This situation is 
strikingly different from that for elliptic equations, such as the time-independent 
Schrodinger equation, where a caustic marks the breakdown of the semiclassical 
approximation. Note that the celebrated turning points of the harmonic oscillator 
are not caustics for the time-dependent problem!) 

With this background understanding we can finish treating the trajectory 
gT5l) : 

• A variant of the calculation leading to (|4.16p shows that, as claimed, S = 
for any trajectory linking y — to x — (necessarily in elapsed time n/u>). 

• As for the harmonic-oscillator kernel, the singularity in A of (|4.15j) is 
harmless; the companion factor e lS l n is effectively there. 

• By continuity from the (purely harmonic) case y = 0, when y < but 
small one would expect a caustic to occur somewhere near x = —y, t = 
tt/uj. Therefore, when the trajectory reemerges from the potential, this 
term of the kernel should be multiplied by a Maslov factor —i. 

To verify this last claim, and to make a start on computing the second (single- 
reflection) term in (|4.9j) . we consider a path that starts at y < at time and ends 
at x > at time t. It must cross the time axis at a time t\ , and from the solution 
of the classical equation one finds 

(4.18) u)xt x + y sm(u(t - h)) = 0, 
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Figure 4. (a) The four possible relations between a diagonal line 
and the principal arc of the sine curve, (b) Resulting division of 
the parameter plane, labeled by intersection numbers. 



which can't be solved by elementary functions. The action is 

2 2 

(4.19) S(x, y, t) = JL + sin (2 W (t -ti)). 

By implicit differentiation of (|4.18[) one can find that 

(4-20) A* = IT 7—7 rr • 

2ti x — ycos{uj{t — ti)) 

There will be a caustic if the denominator of (|4.20|) changes sign. (Since that factor 
arises from dt\/dx, its vanishing says that t\ (hence y) can vary without changing 
x (at least to first order).) 

To investigate further it is helpful to introduce dimensionless variables 

(4.21) T = ut, Sl=uti, <0 = T -Q, <tt, p=-->0. 

x 

Then ([4TTgj) is 

(4.22) = psin(O) + £l-T = f(Cl) (0 < Q < w), 

which can be investigated graphically as the intersection of a trigonometric graph 
and a straight line. The number of intersections can be 0, 2, or 1 (Fig|4|). There is 
one solution if T < it (that is, the straight line hits the axis below the sine curve). 
There are no solutions if p < 1 and T > ir, or if p > 1 and T > T r , where 

(4.23) T* = vV - 1 + cos" 1 (- * 

There are two solutions if p > 1 and 7r < T < T» . Finally, there is one solution 
if the straight line is tangent to the sine curve — that is, f'(Cl) and /(f2) are zero 
simultaneously — which happens when p > 1 and T — T» . Furthermore, 

(4.24) = /'(O) =pcos(fi) + l 

is the condition for the vanishing of the denominator of (|4.20[) . 

Now consider a fixed trajectory with a moving endpoint (that is, fix y and 
ti and let x and t vary). When t w ti , /'(f2) is large and positive (/? — > +cx), 
cos(i7) — > 1). Near the exit point, t w *2 , /'(^) is large and negative (/? — > +00, 
cos(J7) — > — 1). Therefore, every trajectory does pass through a solution of (|4.24[) 
somewhere on its retreat from the potential. 

In future work we hope to complete the calculation of AA for trajectories with 
an endpoint inside the potential. That will enable one to study whether the two 
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conditions in Theorem 14.21 are satisfied in spite of the caustic, as they are for the 
Mehler kernel. If so, then one can tackle the second (single-reflection) term in (|4.9j) 
by concatenating a trajectory of the sort just studied with one of the time-reversed 
type. For given (x, y, t) outside the potential, one must integrate over all (q, r) 
inside the potential for which such a trajectory exists. From the taxonomy of paths 
given above, it is clear that as many as four trajectories can exist, so the term 
GsciQ 1 = J TA. is a sum of four terms, each with a domain of integration that is a 
nontrivial subset of the region < q < oo, < r < t. (Negative q do not contribute, 
because we saw earlier f (|4.13l) and (|4.15[1 ) that AA = there.) 

5. Conclusion 

The seemingly elementary model of a "power wall" has run into several rather 
profound mathematical problems that are worthy of mathematicians' attention. 

(1) What phase shifts S(p) correspond to potentials (or even nonlocal dynam- 
ics) qualitatively worthy of being called "soft walls" ? 

(2) What phase shifts do and don't lead to finite oscillatory integrals (|3.25[) . 
p.26|) . (|3.27|) . (|3.28p . etc.? When the convergence is unstable, can nu- 
merical methods nevertheless be applied to such integrals? Can they be 
analytically recast into more rapidly convergent integrals? 

(3) What can one say in general about caustics in the semiclassical solution of 
time-dependent Schrodinger equations? Are they generically as harmless 
as in the Mehler formula, or as harmful as in elliptic problems? 
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